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Abstract 

^ We define a countable family of sequences of points in the unit square: 

the LS-sequences of points a la Halton. They reveal a very strange 
and interesting behaviour, as well as resonance phenomena, for which 
we have not found an explanation, so far. We conclude with three 

^ open problems. 

Keywords Uniform distribution, Discrepancy, Quasi-Monte Carlo methods. 

a 

Mathematics Subject Classification 11K06, 11K31, 11K38, 11K45. 

1 LS'-sequences in [0, 1[ 

In [I] the first author introduced a countable family of LS-sequences of points 
in [0, 1[, whose construction strongly depends on the LS-sequences of parti- 
t-h Hons of [0, 1[ introduced in the same article. Actually, the former are obtained 

by reordering the left endpoints of the intervals of the latter. 

In this paper we make the first step in the direction of their generalization 
.£h to higher dimension, presenting two possible extensions to dimension 2: the 

LS-point sets a la van der Corput, and the LS-sequences a la Halton. 
03 The interest on low discrepancy sequences of points in dimension two or 

higher is motivated by their application in Quasi-Monte Carlo methods (see 

0)- 

In this section we recall the definition of the LS-sequences of partitions 
and points and their main properties. 

The LS-sequences of partitions are obtained as a particular case of the 
p -refinements introduced by the third author in [10J. 
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Definition 1. Consider any non trivial finite partition p of [0, 1[. The p- 

refinement of a partition tc of [0, 1[ (which will be denoted by p7r) is obtained 
by subdividing only the interval(s) of tc having maximal length homotetically 
to p. Denote by p n 7i the p-refinement of p n ~ 1 n, and by {p n 7r} n >! the sequence 
of successive p ^-refinements. 

If p = {[0, a[, [a, 1[} and u = {[0, 1[} is the trivial partition of [0, 1[, 
the sequence of successive p-refinements is actually the splitting procedure 
introduced by Kakutani [BJ. 

In [TD] it has been proved that the sequence {p n u} n >i is uniformly dis- 
tributed, which means that if p n oo = { [z/i ; [ : 1 ^ * — t n }, the sequence 
satisfies 

lim — V X[a,b[(yi n) ) = b-a, 

i=l 

for every pair of real numbers a, b, with < a < b < 1 (see [10] for the first 
formal definition of this concept). 

An LS-sequence of partitions, denoted by {Pi S } n >i, is obtained as suc- 
cessive pis-refinements of u, where p^s is the partition made by L long 
intervals having length 7 followed by S short ones having length r y 2 , and 7 is 
the positive root of the equation L7 + S7 2 = 1. 

It is clear that each partition p\ s has only long intervals having length 
7™ and short ones having length 7 n+1 . The sequence {t n }n>i of the number 
of intervals of {p£ s)n>i satisfies the difference equation t n = Lt n _i + St n _ 2 
with t = 1 and ti = L + S. 

In pp it has been proved that when S < L the LS-sequence of partitions 
{Pl s) has low discrepancy, namely there exists a constant C > depending 
on L and 5* such that t n D{p r l s ) < C for any n. Here D(p'1 s ) denotes the 
discrepancy of p£ 5 defined by 



D(pts) = sup 

0<a<b<l 



In fact, we have more generally 

Theorem 2. If S < L there exist c\,d x > such that c[ < t n -D(p^s) — c i 
/or any n G N. 

If S = L + 1 there exist c 2 ,c' 2 > such that c' 2 logt n < tnD^p^s) < 
C2 log t n for any fiGN. 
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If S > L + 2 there exist 03,03 > such that d z t^ T < t n D(p2s) < 
c 3^n~ T f or an V n GN, where 1 — r = — > 0. 

In [2] the first author (improving [I]) introduced an algorithm which as- 
sociates to each LS -sequence of partitions a sequence (of points), called LS- 
sequence of points and denoted by {CEg}n>i- The sequences {£i,s}n>i can 
be seen in terms of the representation in base L + S of a suitable subsequence 
of natural numbers. We describe briefly this construction (see [2] for further 
details). 

Any positive integer can be written as 

M 
k=0 

where a k {n) G {0,1,2, . . .,L + S-1} for all < k < M and M = [log L+5 nJ 
(here |_ • J denotes the integer part, as usual). 

Let us denote by Nl,s the infinite set of all positive integers n, ordered by 
magnitude, such that, for each k e N, (ak(n) , a k+1 (n)) {L,L + 1, . . . ,L + 
S - 1} x {1, . . . , L + S - 1}. Define on N Li5 the function 

M 
fe=0 

where dk(n) = afc(n) if < ak(n) < L — 1, while a>k(n) = L + 7(0^(71) — L) if 
L < a k {n) < L + S - 1. 

The sequence {4>L,s( n )}n>o defined on Nl,s is the LS -sequence of points. 

The most important property these sequences show is that whenever the 
LS-sequence of points has low discrepancy, the corresponding LS-sequence 
of points obtained by the algorithm described above has low discrepancy too. 
More precisely, in [1] upper bounds for their discrepancy have been given (see 
[7] for the general theory on uniform distribution and discrepancy). 

Theorem 3. If S < L there exists k\ > such that for any N G N we have 

ND^ls,e L ,s,---,€s) < fologtf. 

If S = L + 1 there exists k 2 ,c' 2 > such that for any N G N we have 
4 log N < N D (ei, 5 , e L ,s, C5) < h log 2 N. 
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If S > L + 2 there exists k^, c' 3 > such that for any N G N we 
have c' 3 N^ < N D fastis , ■ ■ ■ < ^3 N^logN, where 1 - r = 

_i^7) >0 . 

log 7 

2 L/S-sequences of points in the unit square 

In this section we see how the LS -sequences of points can be used to produce 
sequences in the unit square. 

One possibility is to imitate the van der Corput idea [9], combining an 
LS-sequence of points with the points (ordered by magnitude) associated to 
the Knapowski partition j [^r, ^[,l<z<iV}of order N . Hammersely [5] 
extended this definition to higher dimension. 

The other possibility is to put on the two coordinates two different LS- 
sequences of points, imitating what Halton did in jl] when he paired two van 
der Corput sequences having different bases. He proved that whenever these 
bases are coprime, the sequence has low discrepancy in the unit square. 

The first idea produces the following 

Definition 4. For each LS'-sequence of points {££s}n>i; the finite set of 
points 

^,s(A0 = {(^,^), n = l,...,Jv} 

is called LS-point set a la van der Corput of order N in the unit square. 

The main result concerning these LS-points sets a la van der Corput is 
given by the following 

Proposition 5. If S < L there exists C\ > such that for any N e N we 
have N D (p Lj s(N)J < CilogiV. 

If S = L + 1 there exists C<i > such that for any N G N we have 
Nd(p l , s (N)^) <C 2 \og 2 N. 

If S > L + 2 there exists k 3 , c' 3 > such that for any N G N we have 
ND(P L , S (N)^ < CsJV^logJV, where l-r = - l -^>0. 
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Proof. Let us fix a rectangle R = [0, a[x [0, b[ and an LS'-sequence of points 
{££ ff}n>i- If we denote by fi,s one of the three upper bounds appearing in 
Theorem |3j a simple calculation gives. 



N 
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Taking the supremum over all the rectangles R in the unit square, the 
theorem is completely proved. □ 

Remark Comparing the above proposition to Theorem [3j we conclude 
that any LS-point set a la van der Corput has low discrepancy when S < L. 

Let us now give the second and more interesting generalization. 

Definition 6. Given two LS'-sequences of points {££ s,}n>i an d {£Jl 2 &}n>i> 
the sequence 

is called LS -sequence of points a la Halton in the unit square. 
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Halton proved that two van der Corput sequences generate a uniformly 
distributed sequence in the square if and only if they have coprime bases. 

The situation is not clear at all if we use pairs of LS -sequences of points, 
as it can be seen in Figures 1 and 2, where numerical experiments suggest 
that there is no uniform distribution, while Figures 3 and 4 suggest that for 
those pairs of LS -sequences of points we do have uniform distribution. 

The following proposition may give an explanation for the "resonance" 
which appears in Figure 1. 
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Proposition 7. If for each n G N we denote by t n the total number of 
intervals of p\ x and by t n the total number of intervals of the partition p\ x , 
then we have that 

t n = t 3n for all n GN . 

Proof. We use induction on n > 1. If n = 1 we have t' x = 5 = t%. 

Assume t n _ 1 = *3( n _i) holds. The relations t n = t n _i + t n _2 and t n = 
4i„_i + t n _ 2 implies t 3n = 4t 3n _ 3 + t 3n ^ 6 . From the inductive assumption 
we have 4t 3 ( n _i) + t 3 ^ 2 ) = ^t n -i + ^n-2> an d therefore the proposition is 
proved. □ 

Open problems 

1. When is {^£ l s 2 }n>i uniformly distributed? 

2. What is its discrepancy? 

3. What happens in higher dimension? 
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